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This research seeks to assess and improve the accuracy of microphysical assump- 
tions used in satellite passive microwave radiative transfer models and retrieval algo- 
rithms by exploiting complementary observations from satellite radiometers, such as 
TRMM/AMSR-E/GPM, and coincident aircraft instruments, such as the next gen- 
eration precipitation radar (PR-2). We focus in particular on aircraft data obtained 
during the Wakasa Bay field experiment, Japan 2003, pertaining to surface snowfall 
events. The observations of vertical profiles of reflectivity and Doppler-derived fall 
speeds are used in conjunction with the radiometric measurements to identify 1-D 
profiles of precipitation particle types, sizes, and concentrations that are consistent 
with the observations. 


Introduction 

In the mid- to high-latitudes, the spatial structure and microphysics of precipitating 
clouds differ significantly from that of tropical and warm-cloud precipitation (1). 
The variety of precipitation particles, low precipitation rates, surface snow-cover, 
and low freezing level heights all contribute to additional uncertainties involving 
passive microwave (PMW) precipitation retrievals (2). To date, very few coincident 
aircraft and microwave observations have been made of mid-to-high latitude cold- 
cloud precipitation. 

The key research questions to be briefly addressed in this paper are: 

1. What are the ranges of microphysical features, such as particle shapes, fall 
speeds, composition, and sizes, of mid- to high-latitude cold-cloud precipita- 
tion? 

2. What are the primary physical relationships between the microphysics and 
microwave radiative transfer? 

3. To what degree can independent observations, for example, aircraft Doppler 
radar, be used to improve the accuracy of PMW precipitation retrievals? 

4. To what extent and under what conditions can existing PMW methods reliably 
infer details of these precipitating systems? 

We employ aircraft observations of snowfall obtained during the Wakasa Bay ex- 
periment using data from the following instruments: 2nd generation Precipitation 
Doppler Radar (PR-2) operating at 13.4 GHz and 35.6 GHz; Polarimetric Scan- 
ning Radiometer (PSR) at 10.7, 18.7, 37.0, and 89.0 GHz; Airborne Cloud Radar 



(ACR) at 94 GHz; and the Millimeter-wave Imaging Radiometer (MIR) at 89, 150, 
183.3±(1,3,7), 220, and 340 GHz. Satellite observations from the Advanced Mi- 
crowave Scanning Radiometer for EOS (AMSR/E) are also used when reasonable 
spatial and temporal co-location is possible. 

Information inferred from the aircraft observations are used to constrain assumptions 
and parameterizations in a forward model (i.e., simulation of microwave radiances), 
which in turn imposes constraints on the set of solutions computed by a precipitation 
retrieval algorithm (i.e., the inverse of forward model). 

Simulation / Forward Model 

The forward model provides a basis for simulating passive microwave brightness 
temperatures (TBs) and radar reflectivities for a 1-D vertical precipitation profile. 
The three primary components of the forward model are the: physical model , hy- 
drometeor model , and radiative transfer model (RTM). 

1-D Physical Model : The parametric, one dimensional physical model described in 
(3) is employed with modifications. A number of reasonable constraints are applied 
to limit the number of free parameters describing the 1-D column of precipitation. 
Figure (la) illustrates an example of precipitation rate as a function of height for 
pristine snowflakes and aggregates (70% ice by volume). Figure (lb) shows the 3 
GHz simulated radar reflectivity and the Sekhon Srivastava Z-R relationship for 
comparison. 

Hydrometeor Model : Precipitation hydrometeors axe most commonly modeled as ho- 
mogeneous dielectric spheres, so that standard Mie codes may be utilized to compute 
local radiative transfer properties, such as single scattering albedo (wq), asymme- 
try parameter ((g)), and mass extinction coefficient ((fc ex t))- See, for example, (3). 
Each particle is simulated as an effectively homogeneous dielectric sphere of mass 
M composed of a blend of liquid, ice, and air; with volume fractions F such that 
Fuq+Fice+Fair = 1. The density is an average of the individual components (neglect- 
ing the mass of air) where p ^ F\[ q p\\ q + Fi ce Pice- The melted-equivalent particle 
diameter is D — (3M/(4p\ [q )) 1/3 . The physical diameter is D f = D (pu q /p) 1/3 . 
Averaging of the dielectric function uses the effective medium approximation, and 
generalized to three components. 

The distribution of particle sizes is the exponential distribution, given by N(D) = 
No exp[— A £)]. Integration over all diameters yields the total number of particles 
per unit volume of air. The total mass of water per unit volume of air is, therefore, 
W = J 0 °° D s N(D)dD. In the absence of vertical motion, the precipitation 

rate is 1Z = W(v)/p\\ q , where (v) is the mass- weighted average of the particles’ 
terminal fallspeed. For dry snowflakes, observed fall speeds are on the order of 
1 m/s. Aggregates and graupel range up to about 3 m/s (cf., Rogers and Yau, 
Figure 9.7). 

1-D Radiative Transfer Model: The radiative transfer model used for this research is 
the plane-parallel, fully polarized, adding-doubling radiative transfer model, RT4 (4). 



It incorporates surface scattering/emission, and multiple independent layers. Exam- 
ple brightness temperature results for the 1-D model are presented in Figure (lc). 
TB values at common microwave channels are plotted for clear sky (open circles, 
water vapor attenuation only) and a precipitating cloud (Figure la). The brightness 
temperature depression at 89V and 150V (solid circle and line) indicates scattering 
by frozen precipitation, wheres lower frequency microwaves are relatively insensitive 
to scattering by frozen hydrometeors. 
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Figure 1: Simulated 1-D profile for snowfall rate, radar reflectivity, and brightness 
temperatures (51°, over ocean) from the forward model. 

Observations 

Figure 2 shows an example snow case observed by the PR-2 over ocean during 
the Wakasa bay field experiment. The 13.4 GHz and 35.6 GHz reflectivities, nadir 
Doppler velocity, and linear depolarization ratio (at 13.4 GHz) are plotted versus 
height. 4-5 distinct snow storms are apparent in this image: no radar brightband is 
visible and Doppler velocities are mostly confined to the 1-3 m/s range - consistent 
with expected values for snowflakes and aggregates. The 13.4 GHz channel is sen- 
sitive to larger particles but is less susceptible to attenuation through clouds and 
water vapor. Smaller particles scatter microwaves from the 35.6 GHz channel, how- 
ever, attenuation is greater and the reflectivities are generally smaller in comparison 
to the 13.4 GHz channel. 

Other instruments, such as the ACR, MIR, and PSR provide coincident and com- 
plementary datasets, and will be incorporated into the research as necessary and 
available. 


Comparison Plan and Discussion 

PR-2 retrievals of particle size distribution (PSD) parameters are performed us- 
ing the methods of (5), and will be incorporated into a database of simulated 1-D 
model profiles, brightness temperatures, and simulated PR-2 radar profiles - an ex- 
ample of which was shown in Figure 1. However, we do not seek to identify specific 
PSD parameters on a case-by-case basis, but rather identify realistic range bounds 
on PSD parameters for “typical” mid-latitude cold-cloud precipitation cases over 






Figure 2: PR-2 Observations - 1/29/2003 starting 03:18:42 UTC 


ocean and land. One can then match database TBs (and associated profiles) with 
satellite-based passive microwave observations to improve the accuracy of precipi- 
tation retrieval algorithms. 
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